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1. Introduction 

In his 1928 Biometrika contribution [30], Wishart introduced the distribution of covari- 
ance matrices of samples from a normally distributed random variable. His contribution 
triggered a lot of research on theory and application of these and other related multivari- 
ate distributions in e.g., multivariate statistics, probability theory and most recently in 
finance and financial mathematics (for an account of the literature, see |1] and the ref- 
erences therein). An important subclass of Wishart distributions arise as pushforwards 
of normal distributions under certain quadratic forms: Let $,i,C,2, . . . be a sequence of 
M'^- valued independent, normally distributed random variables with mean vectors G M"' 
and covariance matrix E. Then 



3 



is Wishart distributed with scale parameter p := k/2, shape parameter a := 2S and 
parameter of non-centrahty u := Yli=i f^if^J ■ ^^e the notation T{p,uj] a) from [T^ for 
the distribution of S which in turn is motivated by Letac and Massam's [16] family of 
Wishart distributions. Note that S is positive semidefinite almost surely, and also u and 
a are positive semidefinite matrices. Considering first the one-dimensional case d = 1 and 
assuming the non-trivial case S 7^ 0, we can calculate the Laplace transform (or quite 
similiarly the characterteristic function) of ^j^J = by mere completion of a square, 



r-2 i / 2 ('I-A'i)^ 

V2^ 



1 /expf-i±|^r,-^V-i.,?fi-^lM, 



Vl + 27rS ^ /27^II^ V 2^ ^ 1 + 2Sm 



g-«(l+2S«)-Vj2 

(1 + 2Sn)i/2 



integral over density of a normal distribution=l 

n > 0. 



Here E[-] denotes the expectation operator on the respective probability space which sup- 
ports the random variables ^j. By the independence of ^j, (j = 1, . . . , fc), we obtain 

E[e-"^] = rrE[e-"«'] = — — , n > 0. 

^ J ii L J l + cru)P 

j=l 

We see that the family of distributions r{p,u; a) is a natural extension of non-central chi- 
square distributions on the one hand, where S = 1, and of gamma distributions on the 
other hand, where a; = 0. Let now d > 1. In the following we denote by the cone of 
positive semidefinite matrices, by tr(74) the trace of a d x d matrix A, and by det(74) the 
matrix determinant. The positive definite matrices are abbreviated as S'^ . 

With this notation, a similar calculation as above (taking into account the non- 
commutativity of the matrix multiplication) yields the formula for the Laplace transform 

det(/ -I- cru) 



Here / denotes the unit d x d matrix, ab denotes the matrix product of matrices a, b and 
is the inverse of a non-degenerate matrix a. 
It is well known that chi-square distributions and gamma distributions exist for all 
shape parameters p > 0. It therefore may be conjectured that the same holds true in 
dimensions d > 2. However, this is not the case. A number of authors from different 
scientific communities (see the references given in [22]) proved that for invertible a, the 
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central Wishart distributions T{p; a) := T{p, u = 0;a) exist if and only if p belongs to the 
Gindikin ensemble, which equals the set 

Arf := {0, 1/2, ...,{d- l)/2} U (^, oo). 

In other words, the right side of eq. f ll.l|) is the Laplace transform of a distribution on 
if and only if p G A^. 

Note that for d > 3 the set consists of a discrete part and a continuous part. Of course, 
r{p = 0; 0") is trivial (the unit mass at the origin). Some authors therefore exclude 
from the Gindikin set. The non-central case u ^ is more complicated, but also more 
interesting. Peddada and Richards [2T] show by using technically complicated but ele- 
mentary arguments involving special functions that if rank(a;) = 1, then p G A^^. The 
general case (arbitrary rank) has been understood completely very recently. While [T6l 
Proposition 2.3] conjectures the same characterization holds for non-central Wishart dis- 
tributions, the author of the present note has shown in [19] that for p < also the 
parameter of non-centrality must be of lower rank, namely rank(a;) < 2p + 1 (see Theorem 
16.11 below). Furthermore, a preliminary version of that paper, [18], conjectured that in 
this case rank(a;) < 2p (subsequently it turned out that the method I use only implies the 
weaker rank condition rank(a;) < 2p+l). Very recently, Letac and Massam [17] confirm 
my conjecture, while they falsify theirs (see Theorem 17. 5p . 

There is a dynamic way to generate noncentral chi-square distributions. Namely, by 
taking a fc-dimensional standard Brownian motion {Bl, . . . ,B^)~^ and some initial value 
y = {yi)i=i ^ ^'^^ that the non-negative stochastic process 

y/r^Bt) is distributed according to T{k/2, x;2'Et), with initial value Xq = x = y^y > 0. 
This follows from the fact that yi + VT,BI are independent, normally distributed random 
variables with mean yi and variance St. Processes constructed this way are termed " 
Square Bessel Processes", and it can be shown that they are well defined also for any 
non- negative parameter p > 0. For S = 1 and 6 = 2p, Pitman and Yor [26] denote this 
class as W^. 

The matrix-variate generalization of these Square Bessel Processes are the so-called 
Wishart processes introduced by Bru P]. Their crucial feature is the affine property: 
Their Laplace transform is exponentially affine in the initial state, Xq = x. A modern way 
of looking at Wishart processes is by considering them as a subclass of affine processes 
on positive semi-definite matrices or subsets thereof, while the traditional way originating 
from Bru and followed by others is of solving certain stochastic differential equations (in 
these notes they will be called Wishart SDKs). These lecture notes try to explain the 
connections between the two viewpoints. See also section [2l 

We also discuss the existence of Lebesgue densities for Wishart distributions as well 
as the existence of transition densities for Wishart processes. Final remarks are on the 
existence of Wishart processes on state-spaces different from the positive semi-definite 
matrices. 
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2. WiSHART SEMIMARTINGALES, WiSHART DISTRIBUTIONS AND WiSHART PROCESSES 

In this section we introduce and comment on the three main objects of this article: 
Wishart semimartingale, their marginal distributions, which are Wishart distribu- 
tions, and the Wishart processes, which in these notes are Markov processes having 
so-called Wishart transition laws. We will show that Wishart semimartingales can be 
realized as solutions to Wishart SDEs, and that Wishart processes are actually Wishart 
semimart ingales . 

Wishart semimartingales. 

Definition 2.1. Let (1], J-", J-^, P) be a filtered probability space satisfying the usual con- 
ditions. Let p>0, pbeadxd matrix and let further a & S^. A continuous semi- 
martingale Xt is called Wishart semimartingale with parameters {a,p,P), if we can write 
Xt = X -\- Dt + Mt, where Xq = x, Dt = /g(2pa + (3Xs + Xs(3~^)ds and Mt is a local 
martingale with quadratic variation 

[Mt^ij, Mt^ki] = / {{Xs)ikaji + {Xs)iiajk + {Xs)jkaii + {Xs)jiaik) ds, (2.1) 
Jo 

It follows immediately that Mt is continuous with Mq = a.s., and Dq = well. 
An important class of Wishart semimartingales are those obtained by certain squares of 
Ornstein-Uhlenbeck processes. These correspond to drift parameters 2p G N: 

Example 2.2. Let p G N/2, m := 2p, a G and P a d x d matrix. Choose a d x d 
matrix Q for which Q^Q = a (there are, in general, arbitrary many choices for Q). For 
i = 1, . . . ,m we define 

Y.,t := e^* (^y, + e-^'Q^dW^^ , t > 0, 

where W = [W^, . . . ,W"^) is a d x m standard Brownian motion, and is an d- 
dimensional column vector of standard Brownian motions, and yi eW^. 

Then Yi is a Gaussian process for every i = 1, . . . ,m. In fact, Yi is an OU-process 
solving the stochastic differential equation 

dYi^t = PY,,t + Q^dWi, F,,o = y^. 

We define the continuous semimartingale Xt := Xli^i ^«,t^jT- Then Xt starts at Xq = 

Yl^i Villi ! have that 

dXt = {2pQ^Q + f3Xt + Xtf3^)dt + dMt, 

where Mt consists of Brownian terms only. A straightforward calculation yields that Mt 
has quadratic variation fl2.ip . where we have to define a = Q^Q. Hence Xt is a Wishart 
semimartingale with parameters {a,p,P). 

Example 2.3. The following is a particular case of the preceding example (Q = I, f3 = 0), 
but written in matrix form. Let W be a d x n matrix valued Brownian motion, where 
n > d. That is, the entries ofW consist of dxn independent standard Brownian motions. 
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Let further x = yy ■ Then, as can be calculated using Ito-calculus, the process Xt : = 
{y + W){y + W)^ is an -valued Wishart semimartingale with 2p = n, a = I,l3 = 0, 
starting at Xq = x. 

The following note aims at those readers, who are already accustomed to Wishart pro- 
cesses in the sense of Bru: 

Remark 2.4. It is not so trivial to write Xt as solution of a Wishart SDEs, which are 
defined below in equation eq. fl2.2p . The main technical problem is to derive from W and Y 
a new Brownian motion B, for which X satisfies the stochastic differential equation fl2.2p . 
For 2p > d + 1, Pfaffel [251 Theorem 4.19] succeeds by using Levy's characterization of 
Brownian motion. For 2p < d+1 one can show this by an appropriate enlargement of the 
underlying probability space. See statement and proof of Lemma 12. 5[ This technical prob- 
lem supports our decision to introduce the notion of Wishart semimartingale through these 
notes. A further and independent motivation is coming from the recent affine processes 
literature, where the notion of affine semimartingale appears [15]. In our case, the affine 
character of Wishart semimartingales is reflected by the instantaneous drift dDt/dt and 
the instantaneous quadratic variation dMt/dt, which are both affine function in the state 
Xf. The second and important affine character of Wishart semimartingales is constituted 
by their exponentially affine Laplace transform, see Lemma 12.91 below. 

Wishart semimartingales are solutions to Wishart SDEs. We now relate this class 
of semimartingales to solutions of certain stochastic differential equations (SDEs). 

Let y/A denote the unique square root of A & . Let Q, /? be real valued dx d matrices 
and p> 0. As Wishart SDE we define the stochastic differential equation 

dXt = ^/XtdBtQ + Q^dB]^/Xt + {2pQ^Q + pXt + XtP^)dt, Xq = x e S^, (2.2) 

where 5 is a standard d x d Brownian motion. 

We understand any solution of (12. 2 p as weak solution, which means that for given Q, 
(3 and p, there exists a filtered probability space (^2, J-", J-'j, P) which supports a pair of 
Ft adapted processes (X, _B), which satisfy the Ito-integral equation (12. 2 p (which, as is 
common, is written in differentials). However, if the probability space as well as B are 
given in advance, then any solution of (12. 2 p is called a strong one. While every strong 
solution yields a weak solution, the converse does not hold in general. An example of 
a stochastic differential equation (SDE) which admits wealfl but not strong solution is 
Tanaka's one-dimensional equation [201 Example 5.3.2] 

dYt = ^gn{Yt)dWt, Yo = 0. 

While the symmetrization in equation (12. 2 p is necessary to guarantee a stochastic evo- 
lution on the space of symmetric d x d matrices, the existence of solutions for (12. 2p is far 
from trivial. In fact, it is not quite straightforward to ensure that for positive semi definite 

"'^Lct X = B he & standard Brownian motion, then dWt ~ sgn{Bt)dBt is a new Brownian motion 
in virtue of Levy's characterization of Brownian motion, and we have dXt = sg'n{Xt)dWt and clearly 
Xq = Bq = because B is standard. 
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starting values Xq = x G 5*^ , a local solution exists: If we assume x G S'j", then a solution 
exists at least for an almost surely strictly positive stopping time T^. This is the first 
hitting time of the boundary, for X. The solution is a strong one, and its existence follows 
from the fact that the matrix square root is locally Lipschitz on S'j". However, in general, 
the interval [0,TJ^ is purely stochastic, i.e., we might have mf^^Q{Tx{u) \ T^^u) > 0} = 0. 
This definitely happens when p < and under the premise that Q is non-degenerate, see 
Lemma 13.41 and Lemma 17. 4[ For p > , it has been shown in [23] that = oo almost 
surely. In our terminology this means in particular that for each with p > 

a Wishart semimartingale exists. When p E [^^, and in special cases, particularly 
when P = and Q is invertible, Graczyk and Malecki [12] provide the existence of strong 
solutions with different methods than this note (see also the first part of these lecture notes 
written by P. Graczyk). 

We allow for starting values x E S~^, hence in particular we allow that the process Xf 
starts at the boundary dS^ of which are precisely the positive semidefinite matrices 
with rank strictly smaller than d. As the square root is not Lipschitz on dS^, standard 
existence results do not apply. We will see however in a moment, how to infer weak 
solutions from the existence of Wishart semimartingales: 

Lemma 2.5. Any solution of the Wishart SDE (12. 2p is a Wishart semimartingale. Con- 
versely, suppose Xt is a Wishart semimartingale. Then there exists an enlargement of 
{Q, J-", J-i, P) which supports a d x d Brownian motion and some d x d matrix Q for which 
Q^Q = a such that Xt is a solution of the Wishart SDE (12. 2 p . 

Proof. Clearly, the solution of the Wishart SDE is an Ito-process with instantaneous drift 
b{Xt) = {2pQ^Q + pXt + XtP'^). By definition, Xt is the sum of a local martingale (the 
Brownian terms) and a process of finite variation (the integrated drift) plus initial value, 
Xt = X + Dt + Mt- Furthermore, writing out the Brownian terms of (12. 2p in coordinates 
(and using Einstein's summation convention, where summation is performed over all indices 
which appear twice), we have 



Hence, using the formal rules d[Bt,ah, Bt^cd\ = if (a, 6) ^ {c,d) and d[Bt^ab, Bt^cd] = dt if 
(a, b) = (c, d), we have 

d[Xt,ijXt^ki] = d[Mt,ij, Mt^ki] = {{Xt)ikaji + {Xt)iiajk + {Xt)jkaii + {Xt)jiaik) dt, 

where we have set a = Q^Q. Hence we see that Xt is a Wishart semimartingale. 

The converse direction is proved in full generality in [1]. To avoid technicalities (which 
only arise in view of the multivariate character of the problem), and to see the essence of 
the problem, we just consider the case d = 1 here. This is also in some way a prelude 
foreplay for what is demonstrated in more generality in section [31 

We have Xt = Xq + Dt + Mt, where Xq = x, dDt = {h + f3Xt)dt, 6 > and d[M, M]t = 
a^Xtdt. 
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Suppose first x > 0, and b > Using Ito-calculus we see that Yt = log(e ^*Xt) 

satisfies 

dVt = -Pdt + —{dXt + dMt) - 77^(y^Xtdt = ^-^dt + dMt/Xt, 

Jit ^^t -^t 

which equals the differentials of a non-negative process plus a continuous local martingale. 
If Xf would hit zero in finite time, then Yt would go to — oo in finite time. Because the 
first summand above is non-negative, this carries over to the second one. But J^X-^dMs 
is actually just a time changed Brownian motion, hence oscillates infinitely often (and a.s.) 
between — oo and +oo. It can not go to — oo in finite time! So we see that Xt is strictly 
positive a.s., and for all t > 0. Now we can invert Xt, and therefore the process Bt defined 
by 

,^ dXt-dDt dXt-{h + pXt)dt dMt ^ ^ 

CLJDt := 1= — = 1= = — 1^=, -Do •= U, 

cry Xg (Jy Xt cry Xt 

is a well defined continuous local martingale and by construction, [Bt,Bt] = t a.s., for 
all t > 0. Levy's continuity theorem applies and yields that Bt is a standard Brownian 
motion. Rewriting the definition of Bt yields that Xt is a solution of the Wishart SDE 

dXt = {b + (3Xt)dt + a^/xldBt, Xq = x. 

In the general case (where Xt may hit zero in finite time, or even start there), one must 
in general enlarge the probability space to obtain Xt as solution of a corresponding Wishart 
SDE. To this end we use the arguments of [271 Theorem V.20.1], which are much simpler in 
the case d = 1. Let (r2,^,^(,P) be an enlargement of the current probability space which 
supports a standard Brownian motion W independent of X. We define the process 

Bt:= [ e,dMs+ [ PsdWs, 
Jo Jo 

where 9 and p are the predictable processes 

■= — 7=lxs>o, Pt ■= lxt=o- 
o-y/Xs 

Then by construction [B, B]t = t and 5 is a continuous local martingale starting at 0. 
Hence, by Levy's characterization, it is standard Brownian motion. Furthermore 

dXt = dMt + dDt = a^/xldBt + (6 + f3Xt)dt, 

which is just the Wishart SDE in the one- dimensional situation. In the multivariate case, 
6 and p are vectors, whose construction is due to [23 Lemma V.20.7]. □ 

One can show with very little effort that M is an martingale, for instance by using 
the fact that Xt is Wishart distributed (see Lemma 12.91) . since the Wishart distribution 
exhibits exponential moments. 



Wishart semimartingales are Wishart distributed. First, we define tlie family of 
Wisliart distributions, wfiich is motivated by tlie derivation of (11. ip : 

Definition 2.6. We define the non-central Wishart distribution r{p,u]a) on the space of 
symmetric d x d matrices -whenever it exists-by its Laplace transform 

£(r(p,u;;a))(M) = (det(J + aM))~^e-*'-("(^+'^")"'"), ueS^, (2.3) 

where p > denotes its shape parameter, a G S*^ is the scale parameter and the parameter 
of non-centrality equals u & . 

Lemma 2.7. Any Wishart distribution r{p,u!;a) is supported on . 

Proof. It suffices to sliow tliat for any v G M'^, we liave tliat tlie push forward 11,, of T{p, u; a) 
under the map 11 : S'd — ?■ M, x i— )■ v'^xv is supported on M_|_. This, in turn, follows from the 
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fact that n*(r(p, (j){d^) is non-centrally gamma distributed: By Proposition 16.2 
may assume a = 21 without loss of generality. In the following we use A as the Laplace 
variable, and we let U be an orthogonal matrix and fi > such that v = /it/ei, where ei 
is the first canonical basis vector of M'^. Accordingly, u' = U~^ujU. The Laplace transform 
of Il^:{T{p,u; a)){dri) equals 

det(/ + 2Ai;i;^)"Pe"*'^(^^^^(^+2^™^)"''^) 

= det(/ + 2A^2eie7)-^e-*'-('^''^^^^(i+2^'^'^^^^)"''^» = (1 + 2A/i2)-2p/Pe-V(n-2AM^)-^-ii^ 

which is the Laplace transform of fi'^X, where X is a non-central chi-square distributed 
random variable with shape parameter 2p and parameter of non-centrality w[i. □ 

Suppose /3 is again a. d x d matrix with real entries, and let a G S^. In the following we 
denote by the flow of the vector field f3x + x/3^ , that is, 

00^ : RxS^ ^S^, (x) := e^*xe^^*. 
Its twofold integral is denoted by 

a^: R+xS^^st af (y) = 2 /* a;f (z/)rfs. 

Using these two functions, we define a curve 0(t, ■) and a matrix- valued curve ip{t,-) 

(j){t, u) = plogdet (^I + naf (a)) , (2.4) 

ij{t, u) = e^^' (u~^ + trf (a)) ~' e"'. (2.5) 
We show now the elementary fact: 
Proposition 2.8. and ip satisfy a system of generalized Riccati equations, namely, 

^{t,u) = 2ptT{a4j{t,u)), 0(O,n) = O, (2.6) 
iP{t,u) = -2uau + ij{t,u)l3 + I3^ij{t,u), ^Ij{0,u) = u. (2.7) 
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Proof. In order to obtain the generalized Riccati equations fl2.6l) - fl2.7p . we differentiate 
the formula fl2.5p for ip by using the fact that for any differentiable matrix-valued curve 
t I—)- a{t) we have ^^"^(t) = — a~^(t)^a(t)a~^(t), see for instance [9], Proposition III. 4. 2 
(ii)]. Formula (12 ■4p is obtained by using the rule | log(det(a(t)) = tT{a-^{t)f^a{t)), see P 
Proposition II. 3. 3 (i)]. □ 

The following is proved in [1], but with different notation, and for solutions to Wishart 
SDEs. The statement, however, is in fact a result concerning the law of Wishart semi- 
martingales: 

Lemma 2.9. Suppose Xt is a Wishart semimartingale with parameters {a,p,l3) starting 
at Xq = X. Then for each t > 0, Xt ~ T{p, cuf (x); af (a)). 

Proof. Let t > and u G . Let {4',ip) be the functions defined by eqs. fl2.4p - fl2.5l) . 
Applying the Ito-formula to the process 

■— Q-4>{t-s,u)-tT:{ip{t-s,u)Xs) < s < t 

and using thereby Proposition 12. 8^ we obtain 

^ = (a,</.(t -s,u)+ tT{dsij{t - s, u)Xs))ds - tr(^(t - s, u)((/3X, + + 2pQ^ Q)ds + dM,)) 

J s 

+ -4 ti^ipit — s, u)aip(t — s, u)Xs) 

= {ds4>{t -s,u) - 2pti{Q^Q^{t - s, u)))ds - tr(?/^(t - s, u)dMs) 

+ ii{Xs{ds^l}{t - s,u) + 2%l){t - s, u)ai){t - s,u) - %l){t - s, + P'^tl^it - s, u))) 

= -tT{ij{t - s,u)dMs), 

where the first two brackets vanish because of equations fl2.6l) - fl2.7p . We conclude that {Js)s 
is a local martingale on [0, t]. Furthermore, since u) > for alH > and ipit, u) G S*^ 
for all t > 0, we have that J is uniformly bounded on [0,t]. Hence J is a true martingale, 
and therefore 

£[6-"^"^* \Xo = x\= E[Jt \Xo = x\ = J(0) = e-'^(*'")-t'^(^(*'")^), 

where we have used that Jt = e"'^''^"^*) (which follows from 0(0, -u) = and ip{0,u) = u). 
The assertion concerning the distribution of Xt now follows from the explicit formulas 
(I2.4p -( P3]) and the very definition of the Wishart distribution (12. 3p . 

For the derivation of the exponentially affine characteristic function on general state- 
spaces, see the proof of [TOl Theorem 2.2], which uses similar arguments. □ 

Wishart processes from the Markovian viewpoint. 

Definition 2.10. A family of distributions {pt{x,dC,) | t > 0,x G S^} which is non- 
centrally Wishart distributed according to 

Ptix,dO = r(p,a;f(x);af(«))(rfO (2.8) 
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is termed Wishart transition function with constant drift parameter p > 0, linear drift 
parameter /3 and diffusion coefficient a & . 

By using the Laplace transform of the Wishart distribution, we obtain that the Laplace 
transform of the laws pt{x, d^) is given by 

/ e-''^^^^pt{x,d0 = (det(/ + af(a)n))"%-*''("^'+'^'("^")""'^"^) (2.9) 
^ ^-^{t,u)-trwt,u)x) ^ ueSa, (2.10) 

where {4>,ijj) are of the same form as in fl2.4p -f l?^ . 
We start with the following observation. 

Lemma 2.11. Any Wishart transition function is a Markovian transition function sup- 
ported on . The associated Markovian semigroup {Pt)t>o defined by 

f^PJ(x):= [_f{OPt{x,dO (2.11) 

is a Feller semigroup, that is, Pt reduces to a strongly continuous contraction semigroup 
acting on Co{S~^), the continuous functions on vanishing at infinity. 

We use the terminology Wishart process for Markov processes with Wishart transition 
function. 

Proof. First, by Lemma 12.71 we know that for all t > 0, x G S^, the laws pt{x,dC,) are 
supported on S^. For any Borel set B, measurability of pt{x,B) in {t,x) holds by con- 
struction (and in view of the continuity of the maps erf (a), (x) in {t,x)) . So the family 
{Pt)t of linear maps defined by (12. lip is well defined on Bb{S^ ), the set of bounded, Borel 
measurable functions on S*^. We only need to show that it gives rise to a semigroup on 

Since the linear hull of the family of exponentials {/u(0 •= exp(— tr(u,^)) | u G S^) is 
dense in the space of continuous functions vanishing at infinity (and therefore ultimately in 
B{S~^), it suffices to show the semigroup property for the exponential functions ^ i— >■ fuiO^ 
u G S'j". Now, since by Proposition 12.81 we have that (0,'?/^) are the unique solutions to 
a system of ordinary differential equations, it follows (from their specific form) that they 
satisfy the so-called semifiow equations 



(pit + s,u) = (pit, u) + 0(s, ipit, u)) 
ip{t + s,u) = ipis, il>{t, u)). 
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Hence we can write 

-<l}{t+s,u)—tT:{ip{t+s,u)x) 



Pt+sfu{x) = _ fu{OPt+s{x,d^) = e 
Jst 



It remains to prove the Feller property. By [281 Proposition III. 2. 4] and using some 
density argument, it suffices to show that 

• Ptfu{.x) G Co(S'j") for all t > 0, and u E S'^ and this can be seen by inspection of 
ilj{t,u), which is strictly positive definite, as well. 

• Ptfu{x) converges pointwise to /^(x) as t — j- 0, which follows immediately from the 
continuity of (f){t, u) and ?/'(t, u) in t. 

□ 

A Markov process with transition laws pt{x,d^) on is called affine if eq. (I2.10p 
holds. Hence it is obvious that 

Lemma 2.12. Wishart processes are affine processes. 
The generator of a Wishart process. 

Lemma 2.13. Let X he a Wishart process on with admissible parameters {p,f3,a). 
Then the associated semigroup {Pt)t>o has infinitesimal genator A acting on C V^A) 
as 

^/(^) = 1 '^'^''^''^ fx^-dx\i + ^'^^^"^ + ""^^ + 2pQ^Q)Vf{x)), (2.12) 

l<i,j,k,l<d ^-^ 

where V/(x) = and Aijki{x) = {xiuaji + xuajk + XjkOtu + Xjiaik). 

There are different possible proofs of this fact. By using the fact that X can be realized as 
solution of a corresponding Wishart SDE Xt starting at Xq = x, one could just determine 
the generator of X by applying the Ito-formula or using general results on Ito-diffusions. 
Another, maybe more elegant way is the following. By the very definitions of the Wishart 
process, we can calculate the pointwise limit 

lim ~ = 

40 t 

^-m^)~trim,u)-u)x) _ I _ tr((^/;(t, u) - u)x) 
= Lfxjlim ; — ; ; -— 

^ ' m -((){t,u) -tT{{lp{t,u) - u)x) t 

= fu{x) {dtcpit, u) - ti{dti!{t, u)x)) \t=o 

= -fu{x) {2ptr{Q~^Qu) - ti{-{2uau + u(3 + (3^u)x)) . (2.13) 
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The convergence actually holds in sup-norm; this essentially follows from the fact that the 
pointwise limit lies in Co{S^) (see [H Proof of Proposition 4.12]). Furthermore, a density 
argument proves that elements of can be suitably approximated by the linear hull of 
exponentials /^(x), m G S'J^. B On the other hand it is readily checked that (12.1 2p evaluated 
at fu{x), u & equals eq. (I2.13p . 

The drift condition. 

Theorem 2.14. Let X be a Wishart process on with parameters [p, /3, a), and suppose 
a 7^ 0. Then we must have p > 

Proof. Any positive Feller semigroup has an infinitesimal generator A which satisfies the 
strong maximum principle. That is, let / G "^(-4.) and f{x) > /(xq) for all x G . Then 
•^f{xo) < 0. (here the following analogy from calculus helps to remember the sign: Let g 
be a twice differentiable function on an interval /CM which has a local maximum at xq. 
Then f"{xo) < 0. If, in addition, xq lies in the interior of /, then /'(xq) = 0. The analogy 
comes from the fact that the generator of a Feller semigroup has a principal symbol which 
is differential operator of second order). In [1] we used the determinant f{x) = det(x) and 
(diagonal) boundary points Xq G dS^ , because / vanishes precisely there. The theory of 
[1] is more general than these notes, so it is enough to use [H Lemma 4.16 and Lemma 
4.17] to prove the assertion. □ 

Remark 2.15. Note that when a = then we have a deterministic motion, because then 
we have that 

^^-triuXt) I _ ^] _ g-tr(«^f(x))^ 

i.e., Xt = ujf{x). From the Wishart SDE point of view, we clearly have 

X{t) =(3X + Xo = X. 

In that case, p can be anything but is superfluous. 

Wishart processes are Wishart semimartingales. So far we did not need to be spe- 
cific about the realization of Wishart processes as stochastic processes; we only looked at 
the Markovian transition function. In order to relate Wishart processes and Wishart semi- 
martingales, we consider for each initial state x, an associated (to the Wishart transition 
function) Markov process X on a filtered probability space (f2, J-", J-^, P^). Such realizations 
exist and are well known. We repeat in the following a little the definitions for Markov 
transition functions and a canonic construction of the associated stochastic process, which 
is then called Markov. In the end of the section we prove that every Wishart process is a 
Wishart semimartingale. 

A (suitably measurable) family of probability laws t H- {pt{x,dC,)) on S^, indexed by 
t > 0, x G S*^ , is called Markovian transition function, if po{x, dC,) = Sx{dC,) (the unit mass 

^In [4] the rapidly decreasing smooth functions S{S^) arc used. For the corresponding Stone- Weierstrass 
Theorem, see [H Theorem B.3] 
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at x) and it satisfies the Chapman-Kolmogorov equations 

Pt+s{x,A) = I ps{^,A)pt{x,dO, s,t>0. (2.14) 



Note that using the function fA{x) = 1a (a^) (the indicator function on the Borel set A), 
we can write 12.141 equivalently in semigroup form 

Pt+jA{x)=Pt{PsfA)){x), 

where the action Pt is defined above in eq. fl2.1ip . We have therefore shown the Chapman- 
Kolmogorov equations for the continuous functions fu in Lemma (2.111 and that's enough 
by some density argument to ensure fl2.14p . 

Now by P, Theorem 1.1], for any initial distribution z/((i^) on there exists a stochastic 
process X on a filtered probabihty space {Q, T , Tt-, P) whose finite-dimensional distributions 
fulfil 

P[X(o) e Ao, X(t0 e Ai, . . . , x(t„) e A„] 



Ao JAi 

This construction is "canonical" in that Q = {3^^^'°°^ (i.e. the space of all possible paths 
with values in S^), the process is just given by the projections onto the t-th coordinate, 
that is 

and the sigma algebra is given by the product sigma algebra 

The filtration is generated by the projections Xf. 

jr^ = a{Xs I < s < t). 

Starting at z/(x) = Sx{dC,), where x G S^, we denote the associated probability measure 
P by P^. Since f2, J-" and J-i where independent of the initial law, we have constructed 
a family of stochastic processes (^2, J-", J-'f, P^) which satisfy the Markov property for all 
bounded Borel measurable functions /, 

E-[/(X,+,) I J-,] = / f{OPt{Xs,dO = E^M/(^*)] 

which holds P^ a.s., and for all t,s > 0. Here we use to denote the expecation operator 
with respect to P^'. 

This is equivalent to the more intuitive statement 

By Lemma 12.111 we also know that X is a Feller process (this is a Markov process with a 
Feller semigroup), which implies in view of [281 Theorem III. 2.7] that X admits a cadlag 
modification. That means for each x G S"^ , we have that the probability law P^ is actually 
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concentrated on the space of paths which are continuous from the right and have hmits 
from the left. Our aim is to show that for each x & S^, the process Xt is a Wishart 
semimartingale. That is continuous, as we know, which will follow a little indirect: 

Proposition 2.16. Let X be a Wishart process. Then for each x, X is a Wishart semi- 
martingale on [Q, J-", J-'f , P^') . 

Proof. Since X is a Feller process, we have by f2E\ Proposition VII. 1.6] that for any 
fuix) = exp(-tr(Mx)) 

:= /„(Xt) - - / Afu{Xs)ds 

Jo 

is an (J-i, P^)-martingale. Hence by [TU Theorem II. 2. 42] we have that X is a {Q, J-", J-t, P)- 
semimartingale, associated to the generator A. The continuity of X follows from the lack of 
a jump-component in the generator (that is the compensator of the jumps of X vanishes). 
As quadratic variation and drift component are evident from the specific form of the 
generator, we see that X is a Wishart semimartingale on J-", J-'^, P). □ 

3. Boundary non-attainment 

Suppose now that Xq = x is positive definite in (12.21) . In view of the standard existence 
and uniqueness result for SDEs-the square root is analytic, hence locally Lipschitz on S^- 
there exists a unique strong solution of the Wishart SDE as long as Xt does not hit the 
boundary. We call this time 

:= inf{t > | det(X) = 0} 

the first hitting time of the boundary. Of course when Tj. = oo, unique strong solutions of 
the Wishart SDE are guaranteed. This is particular the case, when p is large enough. 

Theorem 3.1. Suppose p > Then T^ = oo almost surely. 

This is a special case of [221 Corollary 3.2] but written in the notation of [H chapter 
6]. The motivation for this result has been the introductory work of [1] for Wishart 
processes with Q = I and (3 = (for more detailed comparison with Bru's work, see [231 
Proposition 3.1]). It should be noted that this is a result concerning the support of Wishart 
semimartingales, and the existence of strong solutions is a mere by-product of the latter. 
For affine jump- diffusions on symmetric cones, the corresponding result is [3 Proposition 
6.1]. 

A random time T : — )■ R_|_ is a random variable taking non-negative values. T is called 
a stopping time, if the sets {T < t} are measurable with respect to J^f In our context 
T = Tj. will always be the first hitting time of solutions to Wishart SDEs of the boundary 
dS~^. [0,T) is called stochastic interval. A local martingale Mf on the stochastic interval 
[0,T) is a stochastic process for which there exists an a.s. strictly increasing sequence 
Tn't T such that for each n, the stopped process Mt/s^T„ is an J^^-martingale. 
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MCKean's argument. This result on continuous semimartingales is fundamental for the 
derivation of Theorem 13. 11 To simplify the setting, we shall from now on assume that T > 
almost surely. This assumption actually holds for T^, because any diffusion started in the 
interior of some domain needs a strictly positive time to reach its boundary. 

Lemma 3.2. For a continuous local martingale on [0, T] almost surely either limt-|-T[M(, Mt] 
exists or we have lim sup^-j^y = — lim inf(-|-T = oo. 

One way of obtaining this result is by performing a time-change on A := 
{limti-T[Mf , Mf] = oo} such that M^t becomes a continuous local martingale(on A) with 
quadratic variation t. Then by Levy's characterization of Brownian motion Mj-^ is a Brow- 
nian motion on A, hence we just need to use the pathwise properties of Brownian motion- 
that a.s. oscillates infinitely often between — oo and oo, as t — oo. The appropriate time 
change is Tf := inf {s > | [M„ M,] > t}. 

A stripped-down version of MCKean's argument is the following. A more general for- 
mulation may be found in [231 Proposition 4.3]: 

Proposition 3.3. Let Z be a continuous adapted stochastic process on a stochastic interval 
[0,T) such that Zq > a.s., and T := M{t > | = 0}. Suppose /i : R+ \ {0} ^ M 
satisfies the following 

(i) for t < T we have h{Zt) = /i(^o) + + Pt, where P is a non-negative process and 
M is a continous local martingale on [0,T). 

(ii) lim^^o h{z) = -oo 
Then T = oo a.s. 

Proof. As a consequence of the assumptions h{Zt) J, — oo as t ^ ^- Since P is non-negative, 
we have that Mt -l —oo as t 'I' T. But M is a continuous local martingale on [0, T). In view 
of the preceding lemma this is only possible, when T = oo. □ 

Now we shortly sketch the proof of Theorem 13. 1[ All the details can be found in an old 
(and unpublished) version of the paper [23] on [211 pp. 5-7]. They base on a few more 
Lemmas. 



then after application of Ito's formula [231 Lemma 4.7] and some lines of calculations we 
obtain 



Proof We define for t E [0,T,) 



Zt := det(e-'^ ^e"''*), h{z) 



log(z). 



h{Zt) = h{Zo) + Mt + Pt 



where 




and 
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where is a one- dimensional standard Brownian motion. Hence M is a continuous local 
martingale on [0,Te] and P is non-negative. Proposition 13.31 can be applied and yields 
Tx = oo. □ 



Hitting the boundary. The following shows that Theorem 13.11 does not hold under 
weaker conditions: 

Lemma 3.4. Let f3,Q be d x d matrices, and suppose Q 0. When p < then there 
exists X G such that < oo with positive probability. 

Proof. Assume, for a contradiction, that for all x G S*^ we have that Tj, = oo. By Lemma 
12.51 any solution Xt of the Wishart SDE is a Wishart semimartingale. And by Lemma 
12.91 we have that Xt ~ r(p, cjf (x); cxf (a)) =: pt{x,d^) where a = Q^Q. By definition 
is a Wishart process. By Theorem 12.141 we must have p > a con- 
tradiction. □ 

For a similar and partially stronger result see Lemma 17.41 below. 

In the case that (3 = and <5 = /, [3 Theorem 1.4] asserts that the boundary is hit in 
finite time, when p G (^^, ^2^)- A similar result including general (3 oi Q ^ seems not 
to be known yet. However, we conjecture 

Conjecture 3.5. Let /3, Q be arbitrary dx d matrices, and let p < Any solution of 
the Wishart SDE with initial condition Xq = x G S]( hits the boundary in finite time, that 
is P(T,. < oo) > 0. 

We further conjecture 

Conjecture 3.6. Let (3, Q be arbitrary d x d matrices, and let p < Any solution of 
the Wishart SDE with initial condition Xq = x G Sj" hits the boundary almost surely, that 
IS P(T, < oo) = 1. 

4. Changing the drift 

Let J-", J^ti P) be a filtered probability space satisfying the usual conditions, which 
supports an J-'^-Brownian motion. Girsanov transformations are tools to derive solutions 
to SDEs as follows. We consider for a moment the one-dimensional case. Let Xt be a 
solution of 

dXt = h{Xt)dt + a{Xt)dWt, Xq = x. 

Suppose we actually seek to solve such an equation for an alternative drift 6(-). If a is 
invertible, we can rewrite the above equation as 

dXt = b{Xt)dt+{b{Xt)-b{Xt))dt+a{Xt)dWt = b{Xt)dt+a{Xt) I ^^^'^ 'J^^" dt + dWt 



aiX, 



b(Xt)-b{Xt) 

probability measure Q, the process '~fsds + Wt is a. Brownian motion, then we have 



In the following we abbreviate jt '■= a{Xt) show that under a new 
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achieved our goal. Note the weak character of this solution: The Brownian motion is not 
given in advance, but constructed from the pair {Xt, Wt). 

What we have outlined above is indeed possible; it is a consequence of Girsanov's theo- 
rem, which asserts that if 

ft 
Jo 

is a martingale on [0,T], then Q defined as 

dQ = £{- [ jsdWs)dF 
Jo 

is a probability measure equivalent to P, and Jg* 'jgds + Wt is a Q-Brownian motion on 
[0, T]. The essential problem therefore is to show the martingale property of {Zt)t<T- That 
can be quite tricky. 

Bru PP used the Girsanov theorem to derive solutions of Wishart SDEs with nonzero 
linear drift /3 from SDEs with constant drift only. In special cases she derives solutions 
until the first time the eigenvalues of the process collide. The respective time of collision is 
not dealt with in her work when /3 7^ 0; recent work elaborates on this issue, see Graczyk 
and Malecki [I2]. 

We have already shown the existence of solutions in the preceding chapter when /3 ^ 
under the more stringend condition p > So we do not need the Girsanov theorem 

to create new solutions, and we also never had to care about the collision of eigenvalues. 
But what we can do is to relate solutions with respect to different drift parameters to each 
other: 

Theorem 4.1. Suppose Xt is a solution of a Wishart SDE with parameters {p, P, Q), 
where Q is invertible, and let Xq = x E . For p^ G M and a d x d matrix f3'^ we set 

7i := v^((/3^ - (/3^)^)g^^ + (p-p^)v^^'g^ (4.1) 

and 

Zt:=£(^-J^ tTijtdBt)y (4.2) 

Ifmm{p,p'^) > ^ii, then Zt is a martingale on [0, T], and dSf := ^t + Bt is a Q-Brownian 
motion on [0,T]. Furthermore Xt satisfies the Wishart SDE with parameters {p'^,P'^,Q) 
under Q. 

Remark 4.2. • On the level of Wishart semimartingales, the result translates in a 

statement concerning their absolute continuity, see [3]. 
• The theorem bases on the fact that under the condition mm{p,p'^) > the 

respective Wishart semimartingales do not attain the boundary dS^ in finite time, 
see Theorem 13.11 Note also: it is impossible to define 7t unless Xt G for all 
t < T. 
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• The best known sufficient criterium for £[— J^'jgdBs) to be a martingale is pro- 
vided by Novikov's condition, 

E[eyoM'dt] < oo. 

This condition is hard to check in our context. Furthermore, it fails, in general. 
In fact, for the particular case Q = I , f3 = 0, P'^ = I , p = we have that 

which is infinite for sufficiently large T. To see this we interpret it as the expo- 
nential moment of a new stochastic process (a so-called affine process) {Xt, Yt = 
tr(Xt)) on X M+ whose moment generating equals 

where ip{t) satisfies the ODE 

= mty + \. m = o. 

This is a matrix Riccati differential equation which has explosion in finite time 
(say at t+ > 0), and by the positivity of x we have that tT{xip(t)) 'j" oo as t ^ 
It follows that the moment f l4.3p explodes. 

For more information on the technique of enlargement of the state space and 
calculation of the moment generating function of affine processes, see for instance 
[TU| Proof of Theorem 4.1]. 

Proof. We start with the second part. Under the premise that Z is a true martingale, the 
conclusion of Girsanov's theorem holds and we obtain 

dXt = ^/XtdBtQ + Q^dBj^t + (2pQTg + + XtP'^)dt (4.4) 

= ^/Xt{-itdt + dBt)Q + Q'^i-iJdt + dBj)y/Xt + (2p<^ Q + + Xt(/3«)^)rft 

(4.5) 

= v%rffi?g + gTd(fi^)7v^ + (2p^QTg + /3^^Xi + Xi(/3'«)T)dt (4.6) 

and therefore Xt is a solution of the Wishart SDE on [0, T] with new parameters (p"^, Q) 
under the measure Q. 

It remains to show that Zt is a true martingale. We use the exact arguments as provided 
by the proof of [21 Theorem 1] but adapted to our matrix-valued setting. 

Since p > we have a well defined positive definite solution Xt of the Wishart SDE 
f l4.6p (subject to Xq = x) in view of Theorem 13.11 and therefore the process 'y{Xt) of 
eq. (14. ip is well defined on < t < T. The stochastic exponential Zt given by f l4.2p is a 
strictly positive local martingale, hence it is a supermartingale. To show that it is a true 
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martingale, it suffices to prov^ that 

E[Zt] = 1. (4.7) 
Quite similarly, there also exists a solution Xt of the Wishart SDE 

dXt = \[^tdBtQ + Q'^dBj^/¥t + (2p^ Q^Q + /J^X^ + Xt{^'^f)dt 

subject to the same initial condition Xt = x (note here: we use the desired new drift 
parameters with Q superscripts, but the SDE is driven by the original Brownian motion 
B). This process serves as auxiliary process to show condition (14. 71) . We also can define 
7t(Xt) exactly as in (14. ip . but using Xt instead of Xt. 
We introduce the two sequences of stopping times 

r„ = inf{t>0| 117(^011 >n}AT 

and ^ 

rn = inf{t>0 I MXt)\\ >n}AT. 
These are increasing sequences satisfying 

lim P(r„ = T) = lim P(r„ = T) (4.8) 

n—>oo n—^oo 

because we use the convention that the infimum of an empty set is +oo. For each n > 1 
we define the process 

:= ^{Xt)lt<r„, te[0,T]. (4.9) 
By construction ||7"(X)|p < n'^t, and therefore Novikov's condition 

Ele^o Ki^n'^i^] <exp{n^T/2) 

which let us conclude that 

= £ l\T{Y:dB,?j 

is a martingale and dQ" := Z^dF defines a probability measure equivalent to P for which 
B^ = 'j^ {X)ds + Bt is a d X d standard Brownian motion on [0,T]. Furthermore, for 
each n, the stopped process X" := Xt^^ have the same law under Q" as the stopped 
processes X^ = Xt/\T„ under P. We therefore have 

E[Zt] = hm E[Z^1,„=t] 



= lim E^"[l.„=r] = lim Q"({r„ = T}) 
= lim P({?w = T}) = 1, 

n— >oo 

where the ffist identity follows from monotone convergence (which is applicable because 
the sets r„ = T are increasing in n, and Z^ is a constant sequence along this sequence; 
hence the sequence Z^lr„=T is a monotonically increasing one). □ 



■^Every positive local martingale is a supermartingale and every supermartingale with constant expec- 
tation is a martingale. See [2H1 P-123]. 



21 

5. On the existence of Wishart distributions 

In this section we provide some results concerning the existence of Wishart distributions 
and their densities. To this end, we introduce some further notation. Let a G M and 
k E Nq. The hypergeo metric coefficient {a)k is defined as 





■ ■ ■ [a + k — 1), otherwise 



Let K = {ki, . . . , Kd) G Nq be a multi-index with length \k\ = ki + ■ ■ ■ + k^. The generalized 
((i-dimensional) hypergeometric coefficient {p)^ is given by 



iph = tl(p-lij-l)) 



see for instance [T5| p. 30]. : S'(i — >■ R shall denote the zonal polynomial of order k, 
where k G Nq. There are several equivalent definitions, for instance 6*^(0 equals the nth 
component of (tr(^))*^, (see [131 Definition L5.1]). Hence 

(tr(0)' = E 

\K\=k 

A more abstract definition [21 P- 234] is that 



C^{i):=uj^ K{k-i)dk (5.1) 

Jk£SO{d) 

where dk is the normalized unique Haar measure on the special orthogonal group SO{d), 
and is some normalizing constant. 

Proposition 5.1. Let p G A^^, cr G S*^ and u G . We have: 

(i) //2p G N and ifTa.nk{u) < 2p, then r{p,uj; a) exists. 

(ii) If p > then the right side of fl2.3p is the Laplace transform of a probability 

measure r{p,uj; a) on S*^. 

(iii) In particular, if p > and if a is invertible, then the density ofT{p,u];a) exists 
and is given by 

F(p,w;a)(0 := (det a)-^ e-*''(""'^+"")(det ef""^ (5.2) 



X 




m!(p)«; I Tk{p) 



where we have set a = a{u) := a ^ua ^, q := q{a) = y/a. 
(iv) If a is degenerate, r{p,u]a) is not absolutely continuous with respect to the 
Lebesgue measure on . 
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Proof. Statement is proved by summing up squares of normally distributed M'^-valued ran- 
dom variables, see section [H 

Note that if a G ^j", our definition of non-central Wishart distribution is related to the 
one of [16] in that T{p, u; a) = 'j{p, a~^ua~^; a), the latter being called "general non-central 



Wishart distribution" in [16]. Hence statement (iii) is a consequence of [161 P- 1400]. 
Now for each £ > we regularize a, a and p by setting 

a£:=a + el, := {cr + el)^^uj{a + el)~^, Pe=p + e. 

1- mr\ 

Then for each e > 0, we pick X^, an valued random variable according to [161 Propo- 
sition 2.3] such that 

Letting e — > and using Levy's continuity theorem, we figure that 

(det(/ + au))"^e-*'^("(-^+"")"''^) = lim (det(/ + a^n))"^^ g-tr(n(7+a,«)-i<^) 

£->0 

must be the Laplace transform of some random variable X ~ r(p, cj;(t), to which Xg, 



converges in distribution as e — j- 0. This settles part (ii) Finally, we consider assertion 



(iv) Assume, by contradiction, that r(p, cj; a) has a Lebesgue density, for some a of rank 
r < d. Let X be an S'^-valued random variable distributed according to r{p,u}; a). Since 
linear transformations do not affect the property of having a density and since the non- 
central Wishart family is invariant under linear transformations (this is easy to check), 
we may without loss of generality assume that a = diag(0,/r), where Ir is the r x r unit 
matrix. Consider the projection 

A simple algebraic manipulation yields that the Laplace transform of nr{X) equals 

which is the Laplace transform of the unit mass concentrated at 7Tr{u}). But the pushforward 
of a measure with density under a projection must have a density again. This yields the 
deserved contradiction. □ 

A very important consequence of this statement in combination of the results of section 
[2] is the following existence result 

Corollary 5.2. For all p > a G S*^ and d x d matrices f3, Wishart processes with 
Wishart transition function with parameters {a,p,P) exist. Therefore all Wishart semi- 
martingales with the same parameters, starting at x ^ exist. Similarly, for all Q with 
Q^Q = a and for all x G 5*^ the Wishart SDE admits global weak solutions. 

Proof. Proposition 15.1 15l allows a well defined Wishart transition function of the form (12. 8p . 
This transition function is Markovian by Lemma 12.111 Now we can combine Proposition 
12.161 Lemma 12.51 to obtain the remaining assertions. □ 
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6. A RANK CONDITION FOR NON-CENTRAL WiSHART DISTRIBUTIONS 

Not for all triples {p,u,a) G M+ x x Wishart distributions r{p,u]a) exist. [19] 
shows the following: 

Theorem 6.1. Let (i G N, p > 0, w G S*^. Suppose a G is invertihle. If the right 
side of (12.31) is the Laplace transform of a non-trivial probability measure r(p, cj; a) on , 
then p ^ Ad and rank(aj) < 2p -\- 1. 



This result contradicts the preceding characterization of Letac and Massam [12], where 
no constraint on the non-centrality parameter had been imposed, which we call here rank- 
condition. Motivated by [19], Letac and Massam [17] deliver very recently an even stronger 
result which uses different methods, and fully characterizes the existence and non-existence 
of the non-central Wishart family (see Theorem 17.51 below). 

Theorem 16.11 uses very nicely the construction and properties of Wishart processes, 
but also elementary arguments, such as Levy's continuity theorem. The latter allows to 
conclude, by using the characterization of central Wishart distributions, that p G A^. The 
proof for the rank condition is indirect; we assume, for a contradiction, the existence of a 
single Wishart distribution which violates the rank condition. We then use the exponential 
family of the latter to construct a whole family of Wishart laws, which determine a Wishart 
process on S^. That, in turn, ultimately violates the drift condition of Theorem 12.141 We 
start with a few lemmas. 

Let {p,u,a) G M++ x 5^ ^ such that /x := T{p,u]a) is a probability measure, that 
is, eq. (12. 3p holds. The domain of its moment generating function is defined as 

D{fi) := {ueSdl C^iu) := [_ e-<"'«V(c?0 < oo}, 

which is the maximal domain to which the Laplace transform, originally defined for u G 
only, can be extended. It is well known that D{n) is a convex (hence connected) set, and 
we also know that C D{fi). Clearly {I + au) is invertible if and only if the (symmetric) 
matrix (J + ^u^) is non-degenerate. Using these facts and the defining equation (12. 3p 
we infer that 

D{^) := {n G 5d I (/ + ^u^) G 5+} = -u'^ + 5+, (6.1) 

and therefore -D(/i) is even open. Accordingly, the natural exponential family of /i is the 
family of probability measure^ 

/ exp«)/i(rfO 

We start by stating some key properties of Wishart distributions @: 



F(/i) 



V G + S\ 



^In order to avoid confusions with calculations in the proof of the upcoming proposition, we change 
here from u notation to u, because u denotes the Fourier-Laplace variable in this paper. 

^Some related properties can be found in Letac and Massam |16j . but in a different notation. Letac and 
Massam use instead of r(p, w; ct) the parameterized family 7(p, a; ct), where a; is replaced by a 
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Proposition 6.2. (i) Letp > 0, u E . Suppose X is an -valued random variable 

distributed according to r{p,uj]I). Let q E and set a := . Then qXq ~ 
T{p, quq; crj^. In particular, T{p, u; I) exists if and only if T{p, quq; a) exists. 

(ii) Let p > 0, a E and u; G S*^ such that fi := r(p, /) is a probability measure. 
For V = — I we have that 



T{p, acua; a). (6.3) 



Conversely, if r{p,auja;a) is a well defined probability measure, so is fi, and 
(16.31) holds. In particular, we have that the exponential family generated by /i 
is a Wishart family and equals 

F(/i) = {r(p, crucr, a) \ a e Sj, a^' - I E D(/i)}. 

(iii) Suppose that T{p, ujq] ao) is a probability measure, forp > and uq, E . Then 
we have: 

(a) r(p, tcjo; ctq) is a probability measure for each t > 0. 

(b) If, in addition, ujq is invertible, then r{p,u;a) is a probability measure for 
each u E , a E . 

Proof. Let E be the corresponding expectation operator. By repeated use of the cychc 
property of the trace and by the product formula for the determinant, we have 
^^-{n,<iXci^ ^ ^^^-{guq,x)^ ^ ^^^^j ^ ^^^yi tT{quq{I + quqY^u)) 

= det(J + au)~^ exp(— ii{uq{I + quq)~^q~^qujq)) 
= det(/ + au)~^ exp(— tr(M(J + au)^^quq)), 



which proves assertion (i) Next we show (ii) We note first, that by (16. ip we have that 



V = a — 1 E D[ji). Hence exponential tilting is admissible. Furthermore, we have 

/ e-<"+'''«>r(p,u;;/)(ciO = dei{I + {u + v)y exp{-ti{{u + v){I + u + vY^u)), (6.4) 

and setting v = cr~^ — I we obtain 

I + u + V = o-~^(I + au). 

Hence the first factor on the right side of eq. (16. 4[) is proportional to det(J + au)^^. It 
remains to show that 

- tr((n + v){I + U + v)~'^uj) = c + tr(n(J + o-n)"Va;cr) (6.5) 



Accordingly (j2.3|) can be written in the form 

£(7(p,a;a))(u) = (det(/ + t7u))^Pe-*'-("(^+"")"'"'^"\ u&S^. (6.2) 

Note that this requires a to be invertible. 

^Expressed in geometric language, we say that the pushfoward of r(p, a;;/) under the map ^ h- !■ q^q 
equals T{p, qojq; <j) 
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for some real constant c, because then the right side of (16 ■4p is proportional to the Laplace 
transform of r(p, auja] a). To this end, we do some elementary algebraic manipulations: 



[u 



v){I 



u 



au)\ 



u 



-OJ 



0- 



(J 



I)uj 
I)iO 
I)iO 



UJ 



aula 



OJ 



+ u) 



UJ 



cr 



UJ 



u 



UJ 



auil + au) auj. 



We set now c := tr((a — /)a;) which is the real number we talked about before. Taking trace 
and performing cyclic permutation inside, we obtain (16. 5p . and therefore the idendity (16. 3 p 
is shown. The assertion concerning the exponential family follows by the very definition of 
the latter. 

We may therefore proceed to 



ni 



and 



Let 



which is proved by repeatedly applying 
r(p, cuo; (To) be a probability measure. Then by (ii) , also Tip^aQ^ujoa^^] I) is one. Let qi 
such that qf = (Ti E S^. We may write T{p, aQ^UJoa^^; I) = T{p, qi^{qiO-Q^ujQaQ^qi)q^^] I), 
and by applying (i) , we obtain the pushforward measure T{p, qiaQ^uoaQ^qi, ai). By (ii) 
we have that T{p,qi^aQ^uJoaQ^qi^; I) is a probability measure as well, and once again by 
we infer that for all a G S* 
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^ , T{p, aqi ^ctq ^wqCq ^Qi ^c", 0") is a probability. We use this 
fact to prove both parts of the assertion. Without loss of generality we assume that a is 
non-degenerate, because in the case a G dS^ we may invoking Levy's continuity theorem^. 
\l\fil and a 



Setting qi 

such that Vg'^CiJoO'o^^'Zr^ 



ctq, we see that 



ma 



holds. For G St we choose q\ G 5*^ 



O UJO 



;iii)b 



□ 



which allows to conclude 
Next, we restate the characterization of the central Wishart laws by using |21j : 
Theorem 6.3. Let d>2, a ^ and p > 0. The following are equivalent: 

(i) det(J + au)~^ is the Laplace transform of a probability measure r{p,u; a) on . 

(ii) p G Ad. 

We are prepared to deliver our proof of Theorem 16. It 
Proof. Let p > such that for some ujq G , a G Sj^ , the right side of (12. 3p is the Laplace 

' ... FTT^ /...X 



111 a 



transform of a non-trivial probability measure T{p,ujo;a). By Proposition 16.2 
have that r{p,iJo/n;a) is a probability measure for each n G N. Letting — > oo and 
invoking Levy's continuity theorem, we obtain that T{p; a) is a probability measure. But 
then by the characterization of central Wishart laws, Theorem l6.3|(ii) , we have that p G A^. 
Let now pq G Ad\[^-^, oo), and let us assume, by contradiction, that there exist {ujq, cr) G 

^ rank(ci;o) > 2po + 1 such that T{pq,ujq; a) is a probability measure. Pick now 



'^Strictly speaking, Levy's eontinuity theorem applies to characteristie functions. However, in the 
Wishart case, the right side of (|2.3|) can even be extended to even the Fourier-Laplace transform with 
ease, and by preserving its functional form. 
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Ui G such that u* := Ui + Uq has rank(ci;*) := rank(ci;i) + rank(a;o) = d, and set 
Pi := ^z£^^l(!i!o) , By construction 2pi = rank(a;i), andpi G Ad\[^^, oo). Hence Proposition 
15.11 15] imphes the existence of a non-central Wishart distribution r(pi, wi, a). Note that 
P* '■= Po + Pi ^ -^d \ [^j^7 oo) and that by convolution 

T{p*,uj*, a) := T{po, ujq, a) -k T{pi, ui, a) 



is a probability measure as well. Since u* is of full rank, we have by Proposition 16.21 (iii)b 
that T{p*,u;a) is a probability measure for all (w, a) G (5'^)^. Hence T{p* ,u;ta) exists 
for all {t,uj,a) G M+ x (sj')^. 

We may now construct a Wishart process by picking some a G S*^ \ {0} and declaring a 
Markovian transition function by setting for each [t, x) G M+ x , pt{x, d^) the probability 
measure given by the Laplace transform 



e 



-<"'«Vt(a;, dO = (det(/ + 2ton))"^* gtr(-«{/+2to«)-i x) ^g_g^ 



(cf. fl2.9p for /3 = 0). Hence X is a Wishart process with constant drift parameter 2p*, 
diffusion coefficient a and zero drift /3 = 0. But 2p* ^ (rf — 1), which contradicts Theorem 
12.141 This shows that we indeed must have rank(a;o) < 2po + 1. □ 

7. Existence of Wishart transition densities 

The aim of this section is to fully characterize the existence of transition densities for 
Wishart processes. That is, we investigate whether the transition laws of Wishart processes 
admit a Lebesgue density. 

Theorem 7.1. Let p > The following are equivalent 

(i) pt{x,d^) has a Lebesgue density Ft^xiO' /^'^ ^^^^ (hence all) t > 0. 

(ii) The d x d"^ matrix 

[Q^ \(3Q^ l-'-lf^'-'Ql (7.1) 

has maximal rank. 

Furthermore, if any of the above conditions are satisfied, then Ft^x is C^{S^), for any t > 0. 

Remark 7.2. • Note that in (17. ip the matrix Q may be replaced by any matrix K 

for which K^K = Q^Q = a. This is obvious from the proof of Proposition 17.31 
below. 

• By an inspection of the (Gaussian) transition law of Ornstein-Uhlenbeck processes 
of the form 

t 



where W is a. (i-dimensional standard Brownian motion, one can infer the well 



known result that (ii) characterizes the existence of Lebesgue densities for Yt. In 
fact, by [731 the covariance matrix of Yf is non-degenerate, for each t > 0, so the 
result holds because F is a Gaussian process. 
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• Condition (ii) is well known in linear control theory, and characterizes the control- 
lability of the linear system 

dtx{t) = I3x{t) + Q^u{t), x(0) = xq. 

That is, let T > 0. Then for each x* G M.'^ there exists a control u such that 
x(T) = X*. For more details, see [23 Chapter 3]. 

The following proposition is a well known ingredient in the characterization of control- 
lability of linear systems. For the sake of completeness and as service for the reader, we 
also prove it here. See, for instance the statements [23 3.1 to 3.4] and their proofs. 

Proposition 7.3. The following are equivalent: 

(i) For one (hence any) t > 0, the matrix erf (a) is positive definite. 

(ii) The d x d"^ matrix fl7.ip has maximal rank. 

Proof. By additivity of the integral, it is clear that if erf (a) is positive definite for some 
t > 0, it is for all s > t > 0. 

By Cayley-Hamilton, for each t > there exist numbers aj{t), j = 1, . . . ,d—l such that 

e*^ = 5^a,(t)A^-. (7.2) 

Hence 

fe^'Q^Qe^^'ds = 5^^7,fc(/3^g^)(/3V) (7.3) 

with 

9jk ■= / aj{s)akis)ds, 
Jo 

which by construction yields a positive semidefinite matrix g := {gjk)jk- 

Proof of (i) ^ (ii) : Since cr^{Q^ Q) is positive definite, we have by using eq. (17. 3p that for 



each z G M^{0} it holds 

5^^,fc(^^/3^g^)(^^/3V)>0. 

But g is positive semidefinite, hence the vector with {z~^ (3^Q~^)j^i must be nonzero. Since 
z was an arbitrary nonzero element of W^, we have proved the rank condition (i 



For the reverse implication, we proceed by an indirect argument. Suppose, there exists 
2; 7^ such that for all t > 0, we have that z'^a^{Q~^Q)z = 0. Due to the positivity of the 
integrand 

Qef^^'z = 

for alH > 0, or equivalently, 

w^Qe^^'z = (7.4) 
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for all w G M'^, t > 0. Since the j-th derivative of e^*Q~^ at t = equals {—iy(3^Q^, we 
have by differentiation of eq. (17. 4p that (3^Q'^ z = for all w and therefore 



cannot 

hold. □ 
Proof of Theorem 17. IL 



Proof. We start with the implication (ii) =^ (i) By Proposition I7.3[ we have that for any 



t > 0, the matrix erf (a) is positive definite. By assumption we have that p > and 
comparing the Laplace transform (12. 9 p of pt{x,d^) with the right side of (12. 3p we realize 



that pt{x,d^) ~ T{p,Uf{x);a^{a)). Hence by Proposition 15.11 (iii) we have that pt{x,d^) 
has a Lebesgue density for each t > 0. 

For the converse direction, we proceed by an indirect argument. Assume, for a con- 
tradiction, that the Kalman matrix (17. ip has rank strictly smaller than d. Then by 



Proposition 17.31 erf (a) is degenerate for some t > 0. But then by Proposition ISTinv) 
Pt{x,dC,) ~ r(p, wf (x); erf (a)) is not absolutely continuous with respect to the Lebesgue 
density. 



So we have shown the equivalence of (i) and (ii) The claim concerning the regularity 
of the densities -Ff,x(0 is an immediate consequence of the second part of Proposition 15.11 

' □ 



m 



Hitting the boundary revisited. As application of this section, we prove the following 



which is stronger to some extent than the assertion of Lemma 13.41 

Lemma 7.4. Let f3,Q be d x d matrices, and suppose that the Kalman matrix (17. ip has 
maximal rank. When p < then for any x G S'j' we have for the solution of the Wishart 
SDE that not only < oo with positive probability but also the stochastic interval [0, T^] 
does not contain a deterministic time interval [0,T], T > 0. 

Proof. Assume, for a contradiction, the existence of x G S'^ for which Tj. > T > 0, where 
T is a positive quantity. An adaption of Lemma 12.51 shows that any solution Xt of the 
Wishart SDE on [0,T] is a Wishart semimartingale on [0,T]. And also from the proof of 
Lemma [2.91 we see that Xt ~ r(p, wf (x); erf (a)) for all t <T, where a = Q^Q. But that 
means that a Wishart distribution exists with non-centrality parameter of full rank, and-in 
view of Proposition 17.31 -also with scale parameter of full rank, but the shape parameter 
satisfies p < This is impossible in view of the subsequent statement. □ 

We cite a special case of Letac and Massam's very recent result [17] on necessary con- 
ditions for the parameters of the Wishart distributions. Translated into our notation it 
reads: 

Theorem 7.5. Suppose a and u are invertible. T{p,u;a) can only exist, if p > 

8. Wishart processes on new state spaces 

In a recent work with Cuchiero, Keller- Ressel and Teichmann [S], the class of afiine 
processes on finite-dimensional symmetric cones H have been completely characterized. 



These are closed convex selfdual cones on which the hnear autoniorphisni group acts transitively 
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Those also contain affine diffusion processes such as the Wishart processes. Symmetric 
cones are classified completely [9], therefore one could try to find SDE realizations as the 
Wishart SDEs (12. 2p on S^. However, only in the case of Hermitian matrices the literature 
provides such realizations. In the latter case we let Wi,W2 be two jointly independent 
d xn Brownian motions {n> d), and y be a complex dxn matrix. Then := {y + Wi + 
iW2){y + Wi- iW2y satisfies 

dXt = ^/XtdBt + dBj ^/Xt + 2pldt, Xo = yf, 

with Bt some d x d complex Brownian matrix, i.e. B = Bi + iB2, where -Bi,i?2 are two 
independent d x d standard Brownian motions. Here c denotes the complex conjugate of 
a complex number c and p = d. Demni |6l chapter 2] dicusses this case, calling these 
processes Laguerre process of integer index. For general drift parameters p > d — 1 see the 
Laguerre processes of [HI chapter 4]. 

[S] delivers for the first time a Wishart process on a non-symmetric cone, namely the 
dual Vinberg cone. This cone is given by the five-dimensional subset K C St 



K 



< u 



( 



\ 



a 


hi 






hi 


Cl 





: u is positive semi-definite > 


b2 





C2) 


> 



and any element x E K can be written as a; = X]^=i yiVl ^ where 



yi 



< ..a ^ 







1/2 



/ \ 

1/2,1 
1/2,2 





2/3 



/ \ 

?/3,l 





yu y yu y y 1/3,3/ 

We give here a slightly different, yet fully equivalent construction. Let {Bi, B2, B^) 



be a 



three dimensional standard Brownian motion. We introduce the process Xt := Y2i=i ^i,tYi[^ 
where for i = 1, 2, 3 we set 

yi,t ■= Vi + Buei 

and ei = (1,0,0)"'' denotes the first canonical basis vector. By Example 12.21 (using /3 = 
0,Q = diag(l, 0, 0) and extending Bi to vector-valued Brownian motions) we know that Xt 
is a Wishart semimartingale on , and by construction Xt is supported on K. Hence by 
the second part of Lemma [275] there exists an enlargement of the original probability space 
which supports a 3 x 3 standard Brownian W motion such that Xt is a weak solution of 
the Wishart SDE 

dXt = i^tdWtQ + Q'^dWi^VXt) + 2pQ^Qdt, Xo = xeK, 

where p = 3/2. 

A full understanding of Wishart processes (leave alone general affine processes) on general 
homogenenous cones is not available at the date this manuscript is printed. 
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A final note might be of interest. By using Example 12.21 one obtains Wishart processes 
with non-convex cone state-space 

Dm '■= {u G I rank(n) < m} 

and these have drift parameter p = m/2, possibly smaller then {d — l)/2, thus violating 
the enigmatic drift condition established in Theorem 12.141 which Wishart processes with 
support on must satisfy. Hence there is no way to extend the so constructed afiine 
processes on Dm to its convex hull 5*^^. This also shows that there are more Wishart 
semimartingales on than those which naturally arise from Wishart processes on . 
But these are supported on the strict submanifolds Dm oi S^, m < d. 
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